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1. DESCRIPTION 


Described in this paper is a computer program that detects outliers in a 
univariate data set. This program, called SOD (statistical outliers 
detection), is capable of detecting as many as 19 outliers in a data set. It 
is written in FORTRAN and can be run in either an interactive or batch mode. 

The SOD software consists of a main program and a subroutine. The main 
program (1) reads the data file, (2) writes the initial section of the report, 
and (3) iterates sequentially for testing the next number of potential 
outliers. The subroutine (1) calculates the test-critical values based on the 
number of potential outliers being tested and compares these to the observed 
values of test statistics and (2) gives the number of observations tested, 
total number of observations considered, mean, standard deviation, extreme 
observation, critical value, and computed test statistics. Also, it prints 
the number of observations declared as outliers and their values. 

The number of potential outliers may be specified by the user or selected by 
the program. Though there is no limitation on the number of observations, it 
is not advisable to use it when there are more than 100 observations. 

It is assumed that the set of observations is from a population which has a 
normal distribution. A significance level of 5 percent is assumed in develop- 
ing the statistical test. If one or more observations do not conform to the 
hypothesis that all observations are from a common population, these observa- 
tions are declared as outliers by the statistical test procedure. Refer to 
appendix A for a brief outline of the test procedure. Details of the 
procedure are described in references 4 and 6. The iterative procedure used 
for specifying the number of potential outliers is discussed in appendix A, 
and a program listing is given in appendix B. 

It should be mentioned that the observations declared as outliers are not 
necessarily "bad" data points, but may be indicative of a nonnormial or 
multimodal distribution; hence, outlying observations should, not necessarily 
be rejected, but must be treated cautiously. 
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2. USER GUIDE 


The SOD program is operative on the PDF Support Processor at NASA/JSC, 
Houston . 

To use the program* the signs on to the computer by simultaneously 
depressing the control key (CHTR) and C key to begin the following dialogue; 
MCR HEL [130,1] (carriage return) 

TASK NAME>S0IL M (carriage return) 

YOUR NAME>User‘s name (carriage return)^ 

MCR>PIP S0D.IN=DATA.DAt 2 (carriage return) 

MCR>RUN SOD (escape) 

The report will be written on the line printer. 
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It may be necessary to simultaneously depress the CNTR and C keys. 


^The data should have been keypunched previously and entered into a file 
(DATA.DAT) with the following format: 


TITLE FOR REPORT 
VALUE 1 
VALUE 2 
VALUE 3 
VALUE 4 


Each number should include a decimal point. 

VALUE N 
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3. EXAMPLE 


The following example illustrates the SOD program. 

MCR>HEL [130.1] 

TASK NAME>SOIL M 

YOUR NAME>HORTON 
MCR>PIP SOD.IN=DANIEL.DAT 
MCR>RUM SODS 
DANIEL (1959) 


0.0 

0.028 

-.056 

-.084 

-.098 

.126 

.168 

.196 

.225 

-.253 

.295 

-.309 

.393 

.407 

.421 

.435 

.463 

-.477 

.547 

.660 

.744 

-.744 

-.758 

-.814 

-.814 

-.898 

1.080 

-1.305 

2.147 

- 2.666 

-3.143 


ENTER NUMBER OF ITERATIONS AS A 2 DIGIT NUMBER OR 
ENTER BLANKS FOR DEFAULT = SQUARE OF NUMBER OF POINTS 


TABLE 1.- COMPUTER OUTPUT FROM SOD 
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. . il.2950 

-Q,3.Q»a- 

0.39JC 

0,9070 

0 , 1 0 

0,9350 

0,9630 

-0,97^0 

0.5970 

0,6600 

0.799() 

•0,7990 

-0.7580 

•0 , 8 I ‘JQ 

-0.6190 

•0.S9S0 

1 . 0^*00 

•1.3050 

2 . 1 970 

• 2 , 6 6 6. 0 

•3.1930 






90PTS c 31 



‘ 



Means « -o. 

1317 





STANDifi'O CrvlATION S 

1 . UOri 1 





NUhBCP 0^ 

TOTil. mijmbFp 


ST 4nj04S'0 


65'J 

UHinCAL lESI 

POTE'lTML 

0UTL1E9S 

0PSE9V47 tQA'S 

mF i 

OtV I i T I'JN 

OHSF.sVi riOM 


V4|,ut5 

6 

31 

•0 ■ 1 3 t 7 

1 ,0'i(l 1 

-3,1930 

3,g1 1 1 

3.3025 

6 

SC 

.03) 9 

0,8935 

-2.6600 

3.1239 

2,8205 

6 

29 

r.0995 

0.6932 

2. 1970 

i , 0 1 1 6 

2,6562 

9 

28 

•0.0)51 

0.5759 

- 1 , 3u5o 

2.2917 

2,5568 

6 

27 

n,03?7 

0.5268 

I ,0600 

1,9882 

2,9936 

6 

26 

•0 , 0 0 76 

0,99J(I 

-0.8980 

1 , a 0 0 3 

2,9906 

5 

3) 

•0,13)7 

1 ,111)0 1 

- 3 , 1 9 3 0 

5.0111 

3,2606 

5 

30 

-0,0319 

0.8935 

-2 , 66O J 

5, 1239 

2,8099 

5 

29 

0,0995 

0,6932 

2.1970 

3.0116 

2.5977 

s 

2« 

-0,0)5) 

0.5759 

-1.3050 

2,2917 

2.5509 

5 

27 

0.0327 

0 , 5 2 e e 

1 .0800 

1 , 986? 

2,97.0a _ 

a 

31 

-0.13)7 

1 . 0 u 0 1 

-3, 1930 

5.0111 

3,20/6 

9 

30 

•0 ,03 1 9 

0.8935 

• 2 , 6 0 6 U 

3. 1239 

2,7629 

9 

29 

o.05‘’5 

0.6932 

2. 1 97u 

5.0119 

2.63Up 

9 

28 

-o.O 1 5i 

0.5754 

• 1 , 3o5o 

2.29 1 7 

2,5ooi 

3 

3! 

•0,1317 

) ,0 0 11) 

-3, 1930 

3.0111 

3 , rib's 

3 

30 

- 0 . 03)9 

0,8935 

-2,6660 

S, 1239 

2,7537 

3 

29 

o.05')5 

0.6932 

2. 1 U7u 

3.0116 

2,5539 


outliers s 3 



POINTS. a -3.193 

-2.6a6 

2.197 
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APPENDIX A 

AN ITERATIVE PROCEDURE EOR DETECTING MULTIPLE OUTLIERS 


APPENDIX A 


AN ITERATIVE PROCEDURE FOR DETECTING MULTIPLE OUTLIERS 
A.l THE STATISTICAL TEST 


Given a set of N observations, suppose one wishes to test the null hypothesis 
of no outliers present in the data set against the alternative hypothesis 
that from one to k. outliers exist, where k is specified in advance. This 
may be done by constructing a sequence of subsets of the data, 

{Ai* ^ 2 * ***» > where Aj is the full set of data and the subset is 

formed by deleting from A^ the observation farthest away from the mean of Aj 
(i = 1, 2, •*•, k - 1). For each subset, the extreme studentized deviate 
(ESD) statistic is defined to be the maximum of the absolute values of the 
studentized residuals. (A studentized residual is the deviation of an 
observation from the sample mean divided by the sample standard deviation.) 

Let t-j be the ESD from the i^^ subset, A.j . Then corresponding to each t^ is a 
critical value such that either (a) t-j < X^ for 1 < i < k, or (b) t^ > X-j 

for 1 < i < h and t^- < x^ for i > h where 1 < h < k. If (a) occurs, the 

hypothesis is accepted that there are no outliersi in case (b), the data are 
declared to have h outliers, with the observations deleted to form as the 
outliers. 

A. 2 CHOICE FOR k 

The critical values x^- in SOD were constructed -for the 5-percent significance 
level by numerical simulation using normally distributed data. These values 
are not independent of k; in fact, for a fixed value of i, they increase 

monotonically with k. As a consequence, even if the basic data are normally 

distributed, the power of the test against a fixed number of outliers 
decreases as k increases so that grossly overspecifying k may result in the 
failure to detect some or even all outliers. On the other hand, if k is 
underspecified, it is more likely that up to k outliers will be detected. 
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but additional ones will not be found since the test assumes no more than k 
outliers are present. In either case, power is lost by using an inappropriate 
value of k. 


If the underlying distribution of the data is nonnormal, especially it is 
multimodal, the test will have a tendency to find many "outliers"; hence, if 
this situation occurs, even for fairly large k, one should be suspicious about 
the distribution of the data and look for underlying mechanisms which might 
have made the data multimodal or highly skewed. 

The following two examples illustrate some of the above points. 


Example 1 : Soil moisture was measured at a depth interval of 5 to 9 cen- 

timeters for a wheat field near Colby, Kansas, on July 18, 1978. The 
following gravemetric measurements of water content in percentage of dry 
weight were obtained from 17 points within the field: 

5.9, 6.4, 5.6, 7.5, 6.7, 4.0, 5.3, 5.5, 5.5, 3.5, 4.6, 10.5, 5.7, 
7..3, 5.2, 9.7, 4.0 


Two observations, 10.5 and 9.7, are suspicious, and one wants to know whether 
or not they could be regarded as outliers. An application of the test 
procedure with k = 2 declares these two observations as outliers; however, 
looking at the data and then choosing k alters the significance level of the 
test by an unpredictable amount. It is interesting to see what would have 
happened had k been chosen (in advance) as 1, 2, 3, or 4. 


Table A-1 provides the computed ESD test statistic, t^- , and the corresponding 
5-percent critical value, (i = 1, 2, •••, k) for different cases. 


TABLE A-1.- ESD STATISTICS AND CRITICAL VALUES 


Subset 

Extreme 

observation 

ESD 

stati Stic 
(t,) 

Critical value 

k = 1 

k = 2 

k = 3 

11 

Al 

lO.o 

2.365 

2.61 

2.74 

2.86 

2.93 

A2 

9.7 

2.549 


2.39 

2.53 

2.57 

A3 

7.5 

1.722 



2.35 

2.44 

A 4 

7.3 

1.728 




2.34 

Number of outliers 

decl ared 

0 

2 

2 

0 


From this table, it can be seen that, if the number of potential outliers were 
specified as either 1 or 4, none of the observations would have been declared 
as an outlier since tj < for all i. In the othei- two cases, t 2 = 2.549 
exceeds both \£^2 ~ ^'2^3 ~ 2.53; hence, the two observations, 9.7 and 

10.5, would be flagged as outliers. Thus, suspected observations may not be 
flagged as outliers by the test when k is under- or over-specified. 

Example 2 : The following soil moisture observations were obtained for a corn 

field near Colby, Kansas, on July 18, 1978, from the top soil layer (0- to 
1-centimeter interval). In this case, gravemetric measurements of water 
content in percentages of dry weight were taken from 35 randomly selected 
points within the field: 

11.5, 3.2, 19.2, 21.6, 5.7, 24.6, 2.1, 3.4, 4.4, 3.7, 4.2, 7.9, 

7.1, 2.6, 3.5, 8.9, 1.8, 2.4, 6.0, 2.8, 29.2, 29.1, 19.6, 1.4, 

4.4, 4.4, 2.9, 4.7, 3.2, 3.8, 2.6, 4.4, 4.6, 4.7, 4.6 


At the top layer, soil moisture can be affected by a number of heterogeneous 
factors; thus, the observations which appear to be outliers may very well be 
legitimate. 



Numerical ordering results in the following data set: 

1.4, 1.8, 2.1, 2.4, 2.6, 2.6, 2.8, 2.9, 3.2, 3.2, 3.4, 3.5, 3.7, 

3.8, 4.2, 4.4, 4.4, 4.4, 4.4, 4.6, 4.6, 4.7, 4.7, 5.7, 6.0, 7.1, 

7.9, 8.9, 11.5, 19.2, 19.6, 21.6, 24.6, 29.1, 29.2 

A quick glance shows that at least one significant gap occurs — between 11.5 

and 19.2. Six observations are greater than or equal to 19.2. When the test 
procedure is applied assuming k = 6, these six observations are declared out- 
liers. Furthermore, its repeated application with k = 2, 3, •••, 10 resulted 
in every additional extreme observation being flagged as an outlier. The 
flagged observations are 7.1, 7.9, 8.9, 11.5, 19.2, 19.6, 21.6, 24.6, 29.1, 
and 29.2. The first four of these observations should not be regarded as 

outliers since water content in this range was found to be quite reasonable 

for places at higher ground on the particular day of measurement. In this 
case, the blind applications of the test leads to the identification of false 
outliers. It must be recognized that these four observations simply cannot be 
Tumped together with the remaining 25 observations and analyzed usi-' data 
analysis techniques based on normal and/or unimodal models. 

Ideally, one should not look at the data before specifying k. Rosner (1975), 
among others, suggested the use of a certain percentage of the number of 
observations to specify k. Barnett and Lewis (1978) proposed that a 
fractional power of N may be used for k. In the author's own work [Chhikara 
(1979)], the rule of k =VN, to the nearest integer, was employed and often 
proved quite satisfactory. Presently, this rule is extended to safeguard 
against errors of undetected outliers as described in the next section. While 
this modified method undoubtedly alters the signifiance level of the test, it 
still provides a useful device for screening data with no prior information. 

If outliers are detected' using this procedure, one should not blindly accept 
such a declaration; instead, this should be taken as a starting point for 
further investigations about the cause of the suspected observations. 
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A .3 AN ITERATIVE METHOD 


Start with k =\/i^ and compute the ESD test statistics. When the statistics t^- 
and the critical values are compared, one of the following cases arises. 

(a) '*'h - ^h k’ ^ ^i ^ k* ^ 

(b) t. < 1 < i < k 

^k ^ H,k 


The test procedure described in section A.l would declare h outliers in case 

(a) , none in case (b), and k outliers in case (c). It is reasonable to assume 
that all potential outliers in the data were detected in case (a); and, hence, 
no further application of the procedure is needed. On the other hand, as we 
have seen, some outliers may have remained undetected in the other two cases; 
i.e., in case (b), the critical values X^^|^ would have been smaller had a 
lesser number been specified for k; whereas, in case (c), additional extreme 
observations might have been declared as outliers had the test been made for 
more than k extreme observations. This suggests that an iterative procedure 
should be used to decrease k in case (b) and increase k in case (c). When k 
is decreased successively by one, a set of smaller critical values is being 
used, thus increasing the power of the test for declaring outliers. In case 

(b) , the decision to stop iterating is made when t^ > X^^ (0 < h < k) which 

results in h outliers being declared, in case (c), In case (c), k is 
increased sucessively so that more extreme observations are tested as 
potential outliers. The test procedure is iterated until 

^k+j I ^k+j,k+j+l 

Vj+1 \+j+l,k+j+l 

with the number of iterations not to exceed a certain preset limit. The 
number of outliers declared is then set equal to (k + j). 
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The critical values estimated by the SOD computer program are based on 
smoothing functions which fit the Monte Carlo-estimated critical values for 
k = 1, 2, 19, and N < 100. Thus, one can apply the test procedure to 

detect as many as 19 outliers in a data set of up to 100 observations. How- 
ever, as pointed out by Chhikara and Feiveson (1980), one would rarely need to 
consider testing for more than a few outliers. Since, at most, 50 percent of 
a set of observations can be thought of as outliers, the SOD program specifies 
k < min (1,19), where I is the largest integer less than or equal to N/2, for 
the upper bound on the number of iterations. 

When using the iterative procedure to obtain k, the significance level of the 
test increases from 5 percent. Any such increase would depend upon the data 
size; however, it is insignificant for large samples since the critical values 
become insensitive to k in such situations. 

A. 4 EXAMPLES 


Returning to Example 1 in section A.2, it is easily seen that the iterative 
test procedure will start by specifying k = 4 under the -JK rule and then will 
consider k = 3, since no outliers are declared in the first case. Because 
^2 ^ ^2 3’ declare two outliers for the data in Example 1. 

•In the case of Example 2, the procedure will continue iterating from the 
initial case of k = 6 to the last case of k = 11, declaring 10 outliers, since 
tio > >>10 11 and t^^ < However, as mentioned earlier, the 10 outlying 

observations are not necessarily outliers and flagging them using the test 
procedure reflects primarily on the data distribution being probably nonnormal 
and at 1 east biomodal . 

Example 3: Daniel (1959) reported the following data consisting of 31 

C 

contrasts in order of absolute value in a 2^ experiment: 

0.000, 0.028, -0.056, -0.084, -0.098, 0.126, 0.168, 0.196, 0.225, 

-0.253, 0.295, -0.309, 0.393, 0.407, 0.421, 0.435, 0.463, -0.477, 

0.547, 0.660, 0.744, -0.744, -0.758, -0.814, -0.814, -0.898, 1.080, 
-1.305, 2.147, -2.666, and -3.143. 





The test procedure started with k = 6 by the >/N rule and declared the last 
three points as outliers. The results, as output from the SOD computer pro- 
gram, are presented in table A-2 and show for each iteration: (a) the number 

of potential outliers specified; (b) the number of observations and the mean 
and standard deviation for each subset, (c) the extreme observation and the 
corresponding computed ESD statistic; (d) the 5-percent critical values; (e) 
the number of outliers declared; and (f) the outliers. 


TABLE A-2.- STATISTICAL OUTLIER DETECTION (EXAMPLE 3) 

Number of 


potential 
outl iers 

Number of 
observations 

Mean 

Standard 

deviation 

Extreme 

observation 

ESD 

Critical 

value 

6 

31 

-0.1317 

1.0001 

-3.1430 

3.0111 

3.3025 


20 

-0.0314 

0.8435 

-2.6660 

3.1234 

2.8205 


29 

0.0595 

0.6932 

2.1470 

3.0116 

2.6562 


28 

-0.0151 

0.5754 

-1 .3050 

2.2417 

2 .5568 


27 

0.0327 

0.5268 

1.0800 

1.9882 

2.4936 


26 

-0.0076 

0.4930 

-0.8980 

1 .8063 

2 .4406 



Number 

of outliers 

= 3 




Outlying 

data points = -3.143, 

-2.666, 2. 

147 


Interestingly, all three 

points, 2 

.147, -2.666, 

and -3.143 [which 

previously 


have been considered highly discordant on the null normal distribution; e.g., 
refer to Barnett (1978)], were flagged here as outliers even though the test- 
ing of a greater number of outliers was considered. Furthermore, the iter- 
ative procedure would have considered smaller values of k, thus increasing 
the power of the test, had these data points not been detected in the first 
1 nstance. 

A .5 CONCLUSION 

The proposed iterative method deals with the problem of specifying the number 
of outliers being tested and should minimize the error that otherwise would 



occur in not detecting outliers when they exist in a data set. Although the 
chances of declaring false outliers may increase, the program still provides 
useful information for data screening prior to subsequent analysis. 

Although the iterative test procedure is considered using the ESD test statis- 
tic, the basic approach can be adapted for any of the test statistics pro- 
posed in the literature; e.g., studentized range (SIR), kurtosis (KUR), and 
R-statistic (RST). Details on these test statistics can be found in Rosner 
( 1975 ). 

It is desirable to determine the actual significance level reached by the 
iterative test procedure, particularly for small sample sizes. Also needed is 
an evaluation of its power againsv the non-null hypothesis of a smaller number 
of outliers. A theoretical solution is probably intractable; therefore, one 
should attempt to make these evaluations by using the Monte Carlo technique. 
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PROGRAM LISTING 


^ORTB*N 

S0B..fXN. 

3004. . 
0002 
0003 

00 04. 
0005 


SOOfc 

0007 

0006 


0009 

3010 

0011 

0012 

3013 

0014 

0015 
3019 


0017 
00 16 


0019 

3020 

0021 

0022 

3023 

0024 

0025 
3026 


0027 

0028 
3029 

0030 

0031 
3032 

0033 

0034 
3035 
0036 

3037 

0036 

3039 
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PROGRAM LISTING 

IV-PLUS V02-SI 
. .„ZTP!SLOCKi/;‘lB .. 

, . LOGICAL*! INFO. 

OIMENstcin TkEtICIo), riTLEtl5) 

CONNON A»G,T,AirK, aK 1 ,a.KK»aK.l,NP,,_H( looo) »0i 1 . 00 Q 1 . I , JnfcUK 1 ,, 
1 EVPPK.f xppki n'STfGT.nurLiE ( luoj 

, . DIMENSION TNf.Q‘(- LOoO 1 . .. . ■ 

DATA TUFTA / ;,4pS , 1.313 , ,3o43 , 1.5302 , ,69.1V , 

1 ,;45268 . ,85333 . ,491 , .2261 , , 1429 / 

C 

C ******** ASSIGN UNITS ******* 

CALL ASSIGN! 1 , I’SOp. JM 1 

CALL ASSIGN(2, iLPi ' ) 

CALL ASSIGf!(3, .’Til ' ) 

C 

C ********* ImITiAlUF constants ******* 

1050 CONTINUE 
S 1 = 0 . 0 
52*0.0 

00 1060 TsI.IOnO 

ocn*o.u 

I)*0 

1080 CONTINUE 
1090 CONTJNHF 

c 

c ************ o^^n OaTa FOOM input fll.*- ********** 

C PEAO title 

PEApi 1 , , 1 SOI f tI aE n ) • 1 = 1 • I '53 
I I So FOOMATf 1544 5 
C 

C «E40 0AT4 POINTS 
1*0 

)?0O continue 

Islti 

PE*On , 1 PSO.ENns I si.^0 )0f 1 5 
1?50 F0PM*TnFl5.55 
si*si+n( I) 
s2*s2Tom*orii 
GO to 1200 
c 
c 

C ♦ CALCUIATF Awn wPItE no, PIS, *'Ean, STanuANu OtvIAIlUN ****** 
1500 continue 
nOPT 5*I-1 
PNsFLOATf NOPTS\ 

4VE*S I /PN 

VAF»(S2 • S1*,ST/Pn) t fPN-1.) 

SDEv=SOPT (V aPI 

‘kRITE(?,155nHTlTLEn)»I = l*l5)>lOCn*I = l»NUPTSJ 
1 550 FORMAT (//, 10 k7i5a 4,//, ' 06 SEP V A T I ONS j ' , ( ^ » 6E 1 1 , 4 j ) 

'■'RITE ( 3, 1 560’ /-TItLF ( n , 1st , 15) 

1560 FORnat cn(,tS4(i,/) 

C 

'‘•RITE 12,1551) nOpTS, »VE, SHtv 
1551 format (/, ' NflPTs s', 15./.' MtANS s'.FlO,4,/, 

1 I standard deviation s' ,F 10.4.///5 
'nRITE 13.15515 NOPTS, 4VF, SO£V 

C 


ORIGINAL PAGE IS 
OP POOR QUALITY 


00<t0 

OOU 

3oa2 

0043 


004U 

0045 

0046 

0047 

0048 

0049 

0050 


0051 

0052 

0053 

0054 

0055 
005h 
0057 


0058 

0059 


'i060 


0061 


0062 

0063 


0064 

0065 

0066 
0067 


0068 


0069 

0070 

0071 


W POOR QDALiry 


C ***** ASSUHf IM1 Tt*L N1JH0EH iTE"ATIONS to Bfc THfc S.0U*Rt KUUT CJ^ 

C OF The NUMBcR of WOIiNTS mITML 1.9 LI *.***. 

T«8CPT{P^) 

NlTtH«T*.S , , . , 

IF (MIT£P .GF, 20i miter * |9 

1700 CONTIMLIJ . _ .. . 

C 

c ***** gtvF user ORjION OF SPEClRtlMG MuHBfcW OP ITtriATIUNJJ ***A*.,. . 
'*RITf.C3, 1710) 

1710 format t' EMTto '^il'^StR OF ITtR*TIU'''S iS A d OICIT NUMsfeK_JH . 

1 ENTER blanks FCIh OEFAULT « SiJUAKF HOOT OF 4U«atR OP HUiNTSt) 

REAOt3. 1 7401 NToIAl - - - 

1740 F0RMAT(I2) 

IFCNTPTiL .GT.O^NlTEPsMTRIiL . . 

1790 WRITE fP, 1 8001 

leoo format (/,< uF total number BTANUAHD 

1 ' FVTRtMf FSo critical test '»/< 

2 ' PpTFmTIAL UBSERVAIIUNS MfeAN UtViATIUN . 

3 ' OHSFRVaTtON VALUES',/,' OUTLIERS',/) 


R£GIm Iterative PROcts.s 
NOITEP s mITER 
KAV I TP s mTTF p 
Tsu.O 

NR ■ MOPTS 

no 190? I = i,mp 

s 0 

call iter ('llTpR, TMf I ^) 


1 900 
1902 

c 

C SHJP llNf aRTER ItFRaTiu-; OATA IS HRImTEp Im SUBROUTiMt 
A R I T F f ? . 1 1 1 1 ) 
nil fcjPmat (Pov) 
f 

C CHECK FUP error !m SllBRUUTiME 

IF (ThFOK .e9. -RoRq^RR.) go TO 9000 


NC further TeStINt. If MUMHER of iTtRATIOMS SPtCIEIEU 
IF (N trial .GT' 0) Gt' TO 80 0 0 


check FriR EKPR ,Gf. theurv nf- last iTEh 

IF ILSTeGT ,F0' NiTfrh ,4^0. maXITH .LT. nLITEH) UJ TO BUlu 

IP (LStfgt ,fg' mi tew ) On rn 20 on 

c 

C EXPP .LT, THEHRY, BaCk i|P 1 ITERAriU-M MmLE-SS MD, of IIEKAIluNS « j, 

1 900 MAXTTP » MAyJ, TR . I 

NITFW e .UAVITR 

IF (NITER ,LE. 0) Go TO flolu 
SO TN l9()n 
C 

C txpp ,Gk, thEHrY nM |As) ttER, IMLWEaSL MiJMbtw OF ITfeRATlO.YS 
2u00 NITER s niTeR ; I 

c 

c REINlTlALTilF Constants and BEGIN ITERATIVE CalCOLATIunS 

NP s nopts 

00 200? I = 1 ,nP 
2002 fK{I) s 0 



(1 


1 



307 £. CALL 1 TE» t.HiTAJ . .. ... . . 

C 

.. C SKIP LINE AFTER ItEI^aTiUN OAT* *'? 1 T Tt.N . u SUttHoy TI Ng. , 

3073 wRITk ( 2 , till) 

C check fop EPPOP Iw Sij^p'uUT INf’ '* 

307a IF (THEOK .Ef?. - 9 pP 99 P 9 t) GO TO 9000 ... .. . . . . .. 

007 5 IF (ElfPP)* .LT.TuEOk , 4 n 0 . E KPPK 1 ,L T , THFOK 1 ) GO TO I 9 i 0 

0076 IF (EXPPK .lT. THeOkj QO TO 2500 . . , 

C 

C &KPP ,GE. THEOpy nN K^H. C.^ECK «tl TH . ITeWATICn ... . .. 

0077 IF (EXPPKI ,LT' ThEokM GO TO 8000 
C 

C EKPP ,CF. TMEOpy fiM kTh, ANO K +1 Th ITERaTIOMS UH 

C EKPPCX) .LE.Th&nPyj-K) Anu E <PP f K ♦ D t GL . THeOPT ( k ♦ U ON Ktl TH irEPATlUN 
3078 2500 IF CNITEPaI .Lf. a| 0 PT.S> 2 ) GO TO 2000 


0079 

C 

C 

c 

GO TO 9200 ‘ . 


3060 

c 

8000 

c 

MTtP z MTF9 ; t 


006 1 
0062 

c 

All 1 0 

9U) 1 

r 

'‘RITE f2.Snit) ^.TTE,P 

FQRh*t ( //, ' nUTLlEHS = ' . 12) 


0063 

oo6a 

0065 

0066 

90 12 
90 13 

IF (NTTfR ,FGi 0) STI.IP 
WRITE f2,90t.3l (OiiTl IF f I) , rs| 

format ( //,' dOPjTs s ',7FUi.j,/,7Pin,i,/,7PM.3j 
STOP 



C 

C 

C AA**A*«««**t EPoil° ^'FSSAGFS *********** 

c 

0067 90 00 “"PITI: rP.POol) ^rTFP 

0089 9U01 FOP'^ATf* FAjal ERpOn ijiiPpjG ITEHATIUM LESi Than ^ N'jN^tPU' 

1 . ' p n pi T s n 

006P STOP 

c 

0090 9200 NOPTS? = NnPT^/p 

0091 i*PlT^ f?, 920t) hirfP. Nf.'PTSE 

0092 920 1 P0P«ATf//,' JTcPAtToI'- mO, ('.ISi'J tKCEL'-JS SO, pOinIS/E 

0093 STOP 
C 

309a E^O 



0001 

3002 __ 

00Q3-_ 


0004 


0005 .... 

0006 
00CI.Z 
0000 
0009 

uolo 
3011 
0012 . 


0013 


0014 

0015 

0016 


0017 

0016 

3019 


3020 


302l 

0022 

0023 

3024 


3025 

0026 

0027 

302P 

0029 

0030 
303 l 

0032 

0033 
3034 
0035 


SUBROUTINE ITEo (NITER, THET*) 

tOi^HON Tj l.» aKl< Li.NP».^ CLQQ.0.1./ .5 U aOiLlil!?! 

1 EXPPK.ExPPKiTt.STEGT, outlie (200) 

01PENS1.0N -ThETaCI), .. . - 

c 

_c t NITI AU 7 r C.DMSTAW.T. /flJL-LiLSI. Ui.aAliaN.AwOE_iXPJ__jSt^ 

LSTE6T * 0 

c calculate constants 

-PN » FLOAT(MP) - - - - 

ARG»PN/2. - FLnATfNi fER) 

, CCNU». '.0CUi532_* FlUa r.(NlT.EH.)...*-J‘lU.AJ(NlJER) 

CONl« ?. - EXoCCONl J 

AKlpTHETA(l) - T.A (2V*EXP t-T^tTA (3) ) t 
1 TNfeTA(S)*LOi;;'FLGAT (MTTfcR) ) 

AKKaTHETAH) • Tri£T A C2) *E)(P (-THE T A ( 3) »P1.0AT (Nl T tP }) 
BKl«ThETAM) • THFTAHO)*LO'G(FLUArCNlTERn 
flK4»THETA{«) - THEtAC91AL0G(FLO4TLNnEPn*C0Nl 

0 > 

w 

c **T**TT«* REGTN iterations ****•»*♦* 

00 3000 IsI.NItF-R 
C 

C SET QI value 

IP f I ,E0. 1 )OI*0' 0 

IF (I.GT, 1 .AmH.t.I T.NlTFP)GI*‘t*P(-fcXPf THETA(6)-TiElAC7) 

1 A FLPAT(I) ) ) 

IP ( I .c<3.NlTEP)nI=l .0 

r 

C SET AKI value 

A4 laOloikii * ( 7 ,-qI ) *AH 1 

IP n .EO. t 1 AkIsjk 1 
IP'n.EP.NlTF®)AKIsAKh 

c 

C SET T VALUE 

IP (NITER. nE.1)t= fPLUArcn-t.) /(PLUATCNlThfi)-l ,) 

C 

C SET HKI VALUE 

TssnwTf n 

Bk Ist*BkK a f t ' -T ) *P k L 
IP n.En.l)RKl=QKI 
IP (I.EQ.N[TeW)bKIsHkk 

c 

C FIND SUN np vAL'ips and values SDUAkeu,. count mjnZERU valuES 
SUNJ«0.0 
SUN2*o.n 

1 COUNT s |i 

DO 2?on .1=1, NP 
IP (NCJ 1 .r.r ,())r,n in 
TCOUNT S ICoUnt * 1 
SUN|=siiNl fO(J) 

SUN2sSIIH?tOM) *Df j) 

??oo continue 

iFricoiiNT.i.T,?Aon To sioo 

COUNTsPLOA'f ICnUNTl 
C 

C COF^PUTfc NFAN, STAmDAmO UEVIATION 


0036 

0037 
0036 


0039. 

0000 

0001 
0002 
0003 
OOOU 
OOOS 


0006 

0007 

0006 

00A9 

0050 


0051 

0052 

0053 
0050 


AV6 t /COL'Nt 

YABl#tSU><2-SUHT*Slj''* ( /CGUNIl /. (oOUNJ-l) 

OEVl«SfJRT(V*B| ^ 

C 

C 7 INO f XTPrMF ViLlip 
BIGbu.O 
00 2000 

^FfHfji^ST.oJCn To 2^00 
OlFFstiPSfnC J)-»VF t ) 

IF(OIFF.GT.RlG>LPC*Tt*J 
IF <OIFF.Gr.PlG>RlG»OlFF 
2000 CONTINIIF 
C 

C SAVF EXPFRIHF^TAL ANQ TheOHICAL VAuUtS 
JF(DtV1 .IT.) .F:2n^FXhPsO.Q 
IFCOfcvl ,Gt. 1 .t:20ifexPpa‘3IG/r)fcVl 

IF ( APg.gF . t ,F -501 A«GLfir. = LOG( APGJ 
IF f AHG.I T.l,E.-50lABGtOG»U.Q 
THfO«AkJ + 9kI*APf;LPG 
C 

C SAVE LOCATIOM nF nliTi lEh 
»-(LOCATElsr 
OUTLTEni = D(| OCiTf j 

“'RITtt?.26''0) MlTfef'’, iCGl'NTfAVtl ,ntvi, OILUCaTEJ, tXPP, fHt.U 
2600 FOP^'ATM^^ . T** Ii»'^''» !'j> 5F ) 1 . 0 ) 


C 



X. 

C SAVE EXpP THtOWV 


[F last Op -.'Ext to last ITtHAliUN * 

0055 

IF 

(I .EG. MTpP) fcXPPx 1 

s 

t»pp 

0056 

IF 

(I .EG. \-lTpP) IHFO-l 


THtn 

0057 

IF 

n .FO. s-ITpO-il txPPK 


£ XHP 

0o5fl 

IF 

( I .£». F!l IpF-l 1 fHF.UK 


THto 

0059 

IF 

(FXPP .GF:. THferl LSItGT 


I 


0060 5000 COMIMif 

0061 «£Tllp^ 

(1062 3100 CO^ T IMI£ 

0q 63 T>-Eljk *-9Q<j'5q'J9' 

0060 P£TUW^ 

0o6«j End 


